How proteins squeeze through polymer networks: a Cartesian lattice study 
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In this paper a lattice model for the diffusional transport of particles in the interphase cell nucleus 
is proposed. The dynamical behaviour of single chains on the lattice is investigated and Rouse scaling 
is verified. Dynamical dense networks are created by a combined version of the bond fluctuation 
method and a Metropolis Monte Carlo algorithm. Semidilute behaviour of the dense chain networks 
is shown. By comparing diffusion of particles in a static and a dynamical chain network , wc 
demonstrate that chain diffusion does not alter the diffusion process of small particles. However, we 
prove that a dynamical network facilitates the transport of large particles. By weighting the mean 
square displacement trajectories of particles in the static chain network data from the dynamical 
network can be reconstructed. Additionally, it is shown that subdiffusive behaviour of particles on 
short time scales results from trapping processes in the crowded environment of the chain network. 
In the presented model a protein with 30 nm diameter has an effective diffusion coefficient of 
1.24 ■ 10 -11 m 2 / s in a chromatin fiber network. 

PACS numbers: 



1. INTRODUCTION 

To accomodate the genome within the the cell nucleus, 
DNA is organized into chromatin. 5% — 12% of the in- 
terphase cell nucleus are filled with a dense network of 
chromatin fibers. In recent years, direct visualization of 
chromatin structure and dynamics in live cells has been 
achieved through significant advances in the field of mi- 
croscopy and molecular biology pQ. The growing aware- 
ness of thermal fluctuations and their impact on the chro- 
matin dynamics has modified the formerly common pic- 
ture of an interphase nucleus which contains randomly 
arranged, but static chromatin. 

In [5] specific chromosome sites were tagged in living 
cells of Saccharomyces cerevisiae with green fluorescent 
protein and the motion of interphase chromatin was mea- 
sured at high resolution and in three dimensions. In [3 
chromosome motion was tracked in Drosophila sperma- 
tocyte nuclei by 3D fluorescence microscopy. A highly 
dynamic chromosome organization was found governed 
by two types of motion: a fast, short range component 
over a 1-2 s time scale and a slower component related to 
long range chromosome motion within the nucleus. The 
motion patterns are consistent with random walks. These 
reports and others jl] reveal that the motion of specific 
chromatin regions is Brownian. However, this motion is 
limited to a nuclear subregion, i.e. a given chromatin 
segment is free to move within only a limited subregion 
of the nucleus. 

Diffusive processes in the cell play a key role in keeping 
the organism alive [5 6 : Molecules transported through 
cell membranes, drugs on their way to respective protein 
receptors, and proteins interacting with specific DNA 
sequences constituting all of the biological functions of 
DNA [3 |8] are diffusion controlled reactions. Further- 
more, proteins approaching their specific target sites on 



DNA are transported by diffusion or even facilitated dif- 
fusion. However, the diffusional transport of molecules 
in the interphase cell nucleus is fundamentally different 
from the normal kind of diffusion which a molecule un- 
dergoes in a homogeneous fluid where the mean square 
displacement of a molecule behaves linear in time t, 
(r 2 (t)) = 6Dt with D as the diffusion coefficient. In the 
cell nucleus the motion of proteins is strongly influenced 
by the presence of the "sticky tangle" of chromatin fibers 
due to steric obstruction and thermal fluctuations of 
the fibers. Fluorescence correlation spectroscopy (FCS) 
studies have shown obstructed diffusion of autofluores- 
cent proteins HO] • Obstructed diffusion or subdiffusion 
is characterized by (r 2 (t)) ~ t a with the anomalous diffu- 
sion exponent a < 1. Other FCS measurements indicate 
that most of the nuclear space is accessible to medium 
sized proteins by simple diffusion and that there is no 
preference for interchromosome territory channels |llj . 
For small proteins up to the size of eGFP-tetramers the 
entire intranuclear chromatin network is freely accessible 
[2"rl] . It is still a matter of extensive discussion to what 
extent macromolecular mobility is affected by structural 
components of the nucleus. This question was adressed 
in our recent work |121 I13j where we investigated the 
diffusional transport of nonbinding and binding proteins 
in static polymer chain networks. However, it is known 
that nuclear constituents much larger than the average 
mesh size of the chromatin network can move in a random 
walk fashion [27]. This implies that the fluctuations of 
the network create transient openings for larger objects. 

Therefore, we have developed a theoretical description 
of the diffusion of nonbinding proteins in a dynamical 
polymer chain network. 

Thcroetical descriptions of dynamical polymer chain 
networks have already been given |14] ; however, they did 
not incorporate realistic structures of chromatin fibers. 
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In the following, the chromatin fibers are approximated 
by flexible polymer chains. Proteins diffusing through the 
fiber network are referred to as particles. We investigate 
in detail how the mean square displacement (msd) of dif- 
fusionally transported particles of various sizes depends 
on the three dimensional geometry of static and dynami- 
cal networks of chromatin fibers and their density in the 
cell nucleus. The results from simulations in both net- 
work types are compared. We then prove by statistical 
analysis of the maximum mean square displacement of 
a particle that a dynamical chain network facilitates the 
diffusional transport of large particles. 

The chromatin network in the cell nucleus during 
the interphase is modeled on a three-dimensional lat- 
tice applying a simplified version of the bond fluctuation 
method [15j . the single site model, in combination with 
a Metropolis Monte Carlo procedure [12] . This mini- 
mizes computational time and effort. The bond fluctua- 
tion model (BFM) was introduced as a new effective al- 
gorithm for the dynamics of polymers by Carmesin et al. 
|15j and provides a very effective means to simulate the 
dynamic behaviour of, e.g., dense polymer melts [17lll8j. 
In section 2 the lattice for chain construction and particle 
diffusion is presented. A brief overview of the particles 
of different sizes is given, reflecting the diffusing proteins 
in the cell nucleus. The simulation of chain dynamics of 
one single chain is described. Afterwards, the simulation 
of chain dynamics of dense chain networks, a combina- 
tion of BFM and MC, is briefly explained. In section 3 
the translational diffusion of the chains is observed. It is 
verified that the dense networks investigated here display 
a semidilutc behaviour and exhibit Rouse dynamics. Fi- 
nally, in section 4 the mean square displacements of par- 
ticles of various sizes in static and dynamical networks 
are compared. It is shown that the diffusion of large 
particles is facilitated by a dynamical chain network and 
how data of the dynamical chain network can be recon- 
structed from data of the static chain network. More- 
over, it is demonstrated that an anomalous behaviour of 
the mean square displacement of the particle on short 
times is caused by a trapping process of the particle in 
the crowded environment of polymer chains. Section 5 
summarizes the obtained results. 
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FIG. 1: Mean square displacement (msd) of one chain 
monomer. Data is obtained by averaging over 50 chain con- 
formations. Msd behaviour for short time scales sketched on 
a logarithmic scale. Dashed line: power law fit, anomaly pa- 
rameter a — 0.45. Inset: msd behaviour for long time scales. 
Straight line: linear fit. 



one occupied lattice site (small), a 2 3 cube (medium), 
and a 3 3 cube (large). The movement of the particle is 
modeled by a random walk. Particles are allowed to visit 
only unoccupied lattice sites. If a particle collides with a 
chain, it is reflected to the last visited lattice site. This 
counts as a time-step, even though the move was unsuc- 
cessful. The initial position is sampled randomly among 
nonoccupied lattice sites. Then, a number of random 
trajectories (typically of the order of 10 3 ) are simulated 
to deliver the ensemble averages. 



2.2. Chain dynamics 



2. THE MODEL SYSTEM 

The system is contained in a 100 x 100 x 100 cubic 
lattice. To prevent boundary effects due to the lattice 
walls, periodic boundary conditions are applied for both 
the particle motion and the chain construction. 



2.1. Simulation of particle movement 

In general studies tracer particles with masses of 27 
to 282 kD were used to study diffusional transport [TU] , 
In our model system we use three different particle sizes: 



The chromatin fibers are modeled as chains con- 
sisting of monomers connected by freely jointed seg- 
ments. Monomers are represented by occupied lattice 
sites. No lattice site is allowed to be occupied more 
than once, otherwise there exists no pair interaction 
between monomers. This corresponds to a real chain 
with excluded volume effect in a good (more precisely: 
athermal) solvent. Allowed bonds between two adjacent 
monomers are taken from the set of all component per- 
mutations and sign inversions P of two bond vectors, 
P(l, 0, 0) U P(l, 1, 0), inducing a bond length of 1 or y/2. 
Chain movements are simulated with a combined version 
of the BFM and MC, which is discussed in detail in [T2"] . 
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FIG. 2: Correlation function pr as a function of time t for 
different chain lengths N. Squares: N = 50, circles: N = 100. 



2.2.1. Single chain dynamics 

Although this paper is dealing with the diffusion of 
particles inside a chain network, it is instructive to first 
analyze the dynamics of single chains, in order to obtain 
the time scales which are required to simulate ensem- 
ble averages of the network. Thus, we investigated the 
dynamical behaviour of single chains of different lengths 
N, N G {10,30,50,100}. Initially, a straight line of N 
connected neighbouring lattice sites is placed into the lat- 
tice and relaxed into its equilibrium conformation. Dur- 
ing one time step (sweep) of the MC procedure, the N 
monomers of the chain sequentially undergo a trial-move. 
If this move hits an occupied lattice site or violates the 
bond length restriction, the chain monomer is set back to 
its previous position and the next monomer is addressed. 



2.2.2. Chain network 



The combined BFM and MC procedure is applied to 
an initial conformation of chains folded into cubes. See 
P2] for a detailed description of the cube construction. 
The relaxed conformation then serves as the initial con- 
formation for investigations of the dynamical network. 
Similarly to the case for single chains, all k x N monomers 
of the k chains, each of length N, are touched during a 
single time step (sweep) of the MC procedure. 



3. RESULTS - CHAIN DYNAMICS 

3.1. Single chain 

The center of mass movement ((x c , m ,(t) — £ c .m.(0)) 2 ) 
over time t is investigated for different single chains of 
chain length N, N e {10,30,50}. x c , m .(t) is the center 
of mass position vector of one chain at time t. Results 
are the average of 50 independent runs. The translational 
diffusion coefficient D tra ns of the chain is obtained from 
the mean squared center of mass displacement : 

6tD trans = ((x c . m .(t) - x c . m .(0)) 2 )). (1) 

First it was verified that the translational diffusion coef- 
ficient D trans decreases roughly inversely proportional to 
the chain length N, 

^trans ^ • (2) 

This is to be expected, because our simulations do not 
consider effects of hydrodynamic correlations (immobile 
solvent approximation) which should lead to the observed 
iV -1 (Rouse) scaling [T^] of the diffusion coefficient. It 
is generally assumed that these hydrodynamic effects are 
playing a less important role inside dense polymeric sys- 
tems, because they are screened on rather short length 
scales. Their omission in the present work is therefore 
justified, considering other, more severe approximations 
which come along with our coarse grained model. 

The mean square displacement of a single chain 
monomer is linear in time on long time scales (Fig. [T| 
inset) 

(r 2 (t)) ~t, *-*■«>. (3) 

On short time scales, however, an anomalous behaviour 
of that quantity is found, 

(r 2 (t)) ~ T, t< 300 (4) 

with an anomaly parameter a of 0.45 (Fig. [T|. It can 
be shown that in case of the ideal chain approximation, 
the anomaly parameter of such a monomer should yield 
approximately 0.5 [20] on short time scales. This is a 
result of the restricted movement of the monomer due 
to its bonds with other chain monomers. The observed 
anomaly parameter of a — 0.45 is a little smaller, which 
we attribute to additional excluded volume restrictions of 
our real (= non-ideal) chain model. On long time scales, 
of course, each monomer has to follow the movement of 
the chain center of mass and hence displays the same 
diffusion behavior. 

For the simulation of dynamical dense chain networks 
it is important to know the relaxation times of the chains 
involved in the simulation, because proper ensemble av- 
erages require simulation times much larger than these 
relaxation times. The slowest relaxation mode of an ideal 
linear chain is the relaxation of its end-to-end vector R, 
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also called Rouse time tr. It is obtained through the 
decay of the time autocorrelation function 



PR(t) = 



(R(0) x R(t)) 

W) 



(5) 



To extract the decay time, it is sufficient to assume a 
single mode exponential decay of the form 



p R {t) = exp 



t 



(6) 



In this way, the Rouse time tr of a single chain with 
length TV = 50 was found as tr ps 2 x 10 4 , and for a chain 
of twice the length, N = 100, the result tr ps 8 x 10 4 was 
obtained (Fig. 111). The Rouse time of an ideal chain 
can be shown to scale with the square of the number of 
monomers I2T1. 



TR~N 2 , 



(7) 



and the same behavior is approximately true for real 
chains. Thus, an important rule of thumb for simulating 
systems of different chain lengths is that the relaxation 
time of the system scales with the square of the chain 
length. 



cd 
E 



CD 
O 

■o 

1i 

E 




2 10 4 



inside an athermal solvent, and hence the fully relaxed 
end-to-end distance of a single chain should yield an av- 
erage of R = dN v , where d p=> 1.2 is the average bond 
length and v the universal scaling exponent, v pa 0.588. 
This yields the overlap concentration of the system in 
athermal solvent 21 



76^— 3 jy4 — Zv 



(8) 



Hence, we obtain <f>*(N = 100) ps 0.03 and <j)*(N = 50) ps 
0.06. Volume fractions below this overlap concentration 
deliver dilute solutions (in which single chains remain un- 
affected by other chains), while volume fractions above 
that threshold deliver semi-dilute solutions. The volume 
fractions in our system were between 0.064 and 0.125, 
above the overlap concentration, hence the system is in 
the semi-dilute regime. The iV _1 Rouse scaling of the dif- 
fusion coefficients remains valid under these conditions, 
but the pre-factors are affected by inter-chain interac- 
tions, i.e. the chain diffusion coefficient differs from the 
single chain case discussed above. 

The center of mass movement of a chain inside dense 
networks with different geometric volumes a, a € 
{6.4%, 8%, 10%, 12.5%} and different chain lengths N, 
N 6 {50, 100} was monitored over the time t. The data 
was averaged over 10 simulation runs. The denser the 
chain network, the slower is the movement of the chains 
and the smaller is the translational diffusion coefficient of 
the chains (Fig. [3]). For a chain network of 6.4% volume 
and chains with length N = 50 a translational diffusion 
coefficient of D tra ns — 0.002 is obtained, whereas chains 
of length N — 100 with the same volume yield a transla- 
tional diffusion coefficient of D trans = 0.00098 (Fig. |3j. 
Hence, the Rouse scaling law as for single chains (see 
equation ([2])) is reproduced, though with a different pre- 
f actor, 
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FIG. 3: Center of mass vs. time t for different chain lengths N 
in a dynamical chain network with different volume fractions. 
Squares: N = 50, (white: D tra ns = 0.0021, black: Dtrans = 
0.0017) circles: N = 100 {Dtrans = 0.001). Blank symbols: 
volume fraction of 6.4%. Solid symbols: volume fraction of 
12.5%. Straight line: linear fit. 



3.2. Chain networks 

In order to characterize the chain network, it is re- 
quired to verify in which density regime the system ex- 
ists. Our chain model is consistent with a real chain 



4. PARTICLE DIFFUSION IN DYNAMICAL 
CHAIN NETWORKS 

In the following the occupation volume of the lattice is 
either computed as the geometric volume or the effective 
volume. The geometric volume is defined as the total 
number of occupied lattice sites. The effective volume 
is defined as the space which is blocked to a particle of 
given size. Hence, the effective volume of the chains on 
the lattice depends on the particle size and the chain 
conformation. For a particle consisting of one occupied 
lattice site, the geometric volume of the chain equals its 
effective volume. 
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4.1. Coupled dynamics of particle and chain 
diffusion 



During one time step of the MC procedure, an entire 
sweep over the chains was carried out, as described in 
section [2~2| and additionally a single move of the walker. 
This setup corresponds to the choice of identical diffusion 
coefficients of the walker (protein) and a single (non- 
connected) statistical chain monomer. In other words, 
a chain monomer and a protein were assumed to be of 
roughly the same size. This does not impose any fun- 
damental restrictions, because more general situations 
could be obtained after a corresponding rescaling of the 
diffusion coefficients. In order to increase the efficiency, 
fifty walkers were moving through the network at the 
same time. They were independent of each other, i.e. 
there was no walker- walker interaction, to simulate the 
ensemble average of single walker systems. The initial 
positions of the walkers were randomly sampled inside 
the lattice, avoiding occupied lattice sites. The simula- 
tions were run for tj^c = 10 5 MC steps; after one run 
was completed, a new set of 50 walkers was initialized. 
Ten runs were averaged in this way. Since the simulation 
time of each single run did exceed the Rouse relaxation 
times of the chains t^c > T R., it was assured that each 
run was featuring a statistically independent set of chain 
conformations. The occpuation volume was varied be- 
tween 6.4% and 12.5%. 



4.2. Physical time scales 

Our simulations yield values for the diffusion coeffi- 
cients D of the walker that are scaled with its diffusion 
coefficient inside an empty lattice, -Do- Consequently, if 
the free diffusion coefficient -Do.exp of a protein is known 
experimentally, it may be multiplied with the diffusion 
coefficient of our simulation to yield the approximate dif- 
fusion coefficient of the protein inside the network. 

Alternately, the Einstein-Stokes relation may be ap- 
plied to estimate real time diffusion coefficients. In our 
coarse grained model, one lattice site approximately cor- 
responds to one Kuhn length of the DNA, which amounts 
to b « lOOnm. Let us assume that we know the physical 
time required for the protein to diffuse a distance equal 
to b, and call it tq. This corresponds to a single time 
step in our random walk simulation, while b stands for 
the step-size. Hence our diffusion law may be interpreted 
as 




6D 



(10) 



where r(t) is a real physical distance, t a physical time 
and D the diffusion coefficient found in the simulation. 
In order to find an approximation for To, we assume the 
protein to diffuse freely on the length scale of one Kuhn- 
lcngth (i.e. one lattice site) and employ the Einstein- 



Stokes relation 



where 



6-Des t , 



k B T 



(11) 



(12) 



Here, fee is the Boltzmann constant, T the absolute tem- 
perature, rj the viscosity of the solvent and R the hydro- 
dynamic radius of the protein. As an example, let us con- 
sider a (quite large) protein with a hydrodynamic radius 



of R = 15nm, T = 300K and rj = lQ-^N ■ , 



With 



k B = 1.3 • 10~ 23 ^ we then obtain D ES ~ 1-4 • 10~ n ^. 
This yields a fundamental time unit of To ~ 1.2 • 10 _4 s. 
Hence, a typical simulation run (10 5 time-steps) would 
correspond to roughly 10 seconds of physical time and 
cover a distance of the order of 10/im. 



4.3. Comparison of particle diffusion in static and 
dynamical polymer networks 

Fig. [4] a displays the dependence of the mean square 
displacement of particles on the time t both in static 
(black symbols) and dynamical (white symbols) networks 
with a geometric chain volume of 6.4%. For the smallest 
particle the mean square displacement shows the same 
behaviour independent of the network dynamics. Fig. [4] 
b displays the dependence of the mean square displace- 
ment of particles on the time t both in static (black sym- 
bols) and dynamical (white symbols) networks with a 
geometric chain volume of 12.5%. 

The mean square displacement of the largest particle in 
a dynamical chain network with an occupation volume of 
6.4% (Fig. |4]a, squares) is higher than in a corresponding 
static chain network. The diffusion coefficients and their 
difference concerning both network types are similar to 
those obtained from simulations with the medium-sized 
particle in a chain network with an occupation volume of 
12.5% (see Fig. |4]b, triangles) due to the same effective 
chain volume of & e ff ~ 52% in both systems. 

The mean square displacement of the largest particle 
in a dynamical chain network with an occupation vol- 
ume of 12.5% (Fig. EJb, squares) is on long time scales 
significantly higher than in a corresponding static chain 
network. This does not imply that each single walker 
is moving faster through the dynamic network. Fig. [5] 
shows that around 5% of the particles (24 out of 500) 
did not reach a maximum mean square displacement be- 
yond 500 (lattice units) 2 (black bars), which means that 
they remained trapped within rather small pockets inside 
the network. In a corresponding dynamical network all 
particles reached a maximum mean square displacement 
larger than 500 (lattice units) 2 (white bars), which nat- 
urally results in a larger averaged value for the diffusion 
coefficient. This finding indicates that most of the walk- 
ers did hardly modify their diffusion inside the dynamic 
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network. A minority of walkers, however, which would 
remain trapped in the static case, becomes released in 
the dynamic case and in this way modify the average 
diffusion coefficient. Therefore, the difference in the dif- 
fusion coefficients is a direct consequence of the averaging 
process during the computation of the mean square dis- 
placement curve. 

In order to verify this finding in a quantitative way, 
the trajectories of the particles in the static chain net- 
work were multiplied with a weighting factor deter- 
mined by the ratio of the both bar lengths (dynami- 
cal/nondynamical). In other words: A walker which 
was moving over a short distance got a lower statistical 
weight than a particle which was moving over a longer 
distance. With this approach, we succeeded to recon- 
struct the data of the dynamical chain network from the 
data of the static chain network (see Fig. |l]b, diamonds). 

It is important to mention that the walker movement 
is not only influenced by the fluctuations of the chain 
network but also by the center of mass diffusion of the 
chains. A pure movement of the walker due to the cen- 
ter of mass diffusion of the chains can be neglected: the 
difference of 0.06 in the diffusion coefficients of the static 
and dynamical chain network with 12.5% volume occu- 
pation (Fig. [4] b) for the largest particle is much larger 
than the translational diffusion coefficient of the chains 
Dtrans = 0.0017 (Fig. [3]) indicating that the movement of 
the walker is mainly influenced by the chain fluctuations. 



4.4. Subdiffusion on short time scales 

Simulations on short time scales with the number of 
MC steps t M c < 5000 were carried out to investigate the 
subdiffusive behaviour of the mean square displacement 
of diffusing particles for short times. 

In a static chain network with a chain volume of 12.5% 
the mean square displacement of the largest particle dis- 
played subdiffusive behaviour with an anomaly parame- 
ter a — 0.81 (Fig. [6] a), in a corresponding dynamical 
chain network a = 0.83 was observed, not significantly 
different when considering the statistical fluctuations. In 
Fig. [6j three different subregimes are observed (dashed 
lines): starting with a diffusion coefficient of D = 0.29, 
the particles become slower after almost 400 MC steps. 
After 800 MC steps the diffusion coefficient decreases to 
D = 0.16. One might suspect that the first collision time 
of the particle could be responsible for this decrease in 
the diffusion coefficient: the particle diffuses freely until 
it encounters the chain network. However, as displayed 
in Fig. [6] b, the majority of particles are colliding with 
the network within the first five time steps. Instead, the 
decreasing diffusion coefficient is caused by the increasing 
number of trapped particles: 16.8% of the investigated 
particles (84 out of 500) did not reach a maximum mean 
square displacement larger than 500. 

Thus, subdiffusion on short time scales is not an effect 
of first collision times of the particle with the chain net- 



work. Instead, subdiffusion results from the fact that on 
short time scales particles are trapped in a crowded envi- 
ronment, hence do not contribute to the large distances 
of the mean square displacement on longer time scales. 



5. DISCUSSION 

The present work was intended to shed some light on 
the diffusional transport of proteins in the dynamical 
chromatin network in the cell nucleus. This process is de- 
scribed with a lattice model of the nucleus, the chromatin 
fibers and a random walk of the proteins. As shown in 
earlier work [T2 [T3] the lattice model has - compared 
to a corresponding continuum model - the advantage of 
being roughly two orders of magnitude faster due to the 
finite number of states. 

To create crowded environments of chromatin fibers in 
the cell nucleus, a self-avoiding random walk chain model 
with excluded volume was used, and chain dynamics was 
implemented through the BFM and MC. The findings in 
Section [3] are hardly surprising to those who are familiar 
with polymer dynamics. The simulation of polymeric 
systems however, is not trivial due to the different time 
scales involeved, and one has to verify that one does in 
fact obtain proper ensemble averages. Thus, we verified 
the known scaling laws and typical relaxation times. 

It may be argued that the chain lengths used in this 
work might be too small to represent a chromatin fiber 
inside a cell nucleus. In reality, the latter would usu- 
ally cover ten thousands of Kuhn monomers and create 
a chain network far inside the concentrated regime (and 
not, as in our case, in the semi-dilute regime). More- 
over, such a chromosome would not exhibit any signifi- 
cant center of mass motion because it remains confined 
inside its territory. In reality, these differences are of lit- 
tle relevance here. The walker does, at any stage of the 
simulation, only "see" a small subsection of the chain. 
Its dynamics is therefore influenced by the fluctuations 
of local chain sections, and hence it is irrelevant how far 
the chain does extend beyond the walker's limited hori- 
zon. In the semi-dilute regime, remote sections of the 
same chain are already uncorrelated, so that a single long 
chain might as well be replaced with several independent 
shorter chains. Additionally, as a result of the Rouse- 
scaling, the chain center of mass diffusion was roughly 
two orders of magnitude slower than the walker's diffu- 
sion and hence it did contribute to its motion only on the 
level of a few percent. 

In several systematic simulations the mean square dis- 
placements of diffusing particles of various sizes were 
computed both in static and dynamical polymer net- 
works. By comparing the results of the two network types 
we could show that chain diffusion does not alter the 
diffusion process of small particles. However, the poly- 
mer network dynamics begins to influence the diffusional 
transport process of particles once the effective volume 
of the chain network is reaching 50% (Fig. [4]). Hence, 
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chain dynamics facilitates the transport of large particles. 
It was demonstrated that the difference in the diffusion 
coefficients in both network types is a direct consequence 
of the averaging process during the computation of the 
mean square displacement. Hence we were able to repro- 
duce the results of the dynamical network after statisti- 
cally reweighting the single walkers contributions to the 
averages obtained within the static network. It might be 
interesting to investigate whether such a procedure could 
be formalized, so that the simulation of dynamical net- 
works could be replaced with static network simulations, 
augmented with a corresponding weight function for each 
trajectory. 

In further simulations anomaly on short time scales is 
found for the msd of a particle in a dynamical fiber net- 
work. We could show that this kind of subdiffusion does 
not result from the first collision time of the particle but 
is caused by a relatively high number of trapped particles 
in the crowded environment of the fiber network (Fig. [6| . 

In recent years, chromatin-binding proteins have been 
extensively characterized due to their functional signifi- 
cance and their dynamics in the living cell has been stud- 
ied by in vivo microscopy techniques [23, 24 . However, it 



is not yet clear how proteins bind to the chromatin fibers 
and displace themselves within the dynamical chromatin 
network pi2"5]. 

We believe that using the model of diffusional trans- 
port of proteins in a dynamical fiber network presented 
here, it is feasible to clarify these questions with an addi- 
tional implementation of a binding process and with the 
help of further systematic simulations. This is ongoing 
work. 
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FIG. 5: Maximum mean square displacment of a particle. 
Number of particles: 500. Black bars: static network, white 
bars: dynamical network. 3x3x3 particle (large), chain 
volume: 12.5%. For a better overview only maximum msd 
values less than 15000 lattice units are shown. 




FIG. 4: Mean square displacement vs. time t for differ- 
ent particle sizes and volume fractions both in static and dy- 
namical networks. White symbols: dynamical chain network, 
black symbols: static chain network. Chain length: N = 50. 
Circles: lxlxl particle (small), triangles: 2x2x2 par- 
ticle (medium-sized), squares: 3x3x3 particle (large), dia- 
monds: 3x3x3 particle in a static network with weighting 
factors. Upper panel: 6.4% chain volume. Circles: D sta uc = 

D dynamical 0.8T. Triangles. D s tatic D dynamical 0.72. 

Squares: D sta tic = 0.49, D dy namicai = 0.53. Lower panel: 
12.5% chain volume. Triangles: D sta tic = 0.44, D dy namicai = 
0.48. Squares: D sta tic = 0.15, D dynamical = 0.21. Diamonds: 

D weighting factor — 0.21 
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FIG. 6: t < 5000. Upper panel: Mean square displacement 
vs. time t on short time scales. Static chain volume: 12.5%, 
3x3x3 particle (large). Middle panel: First contact time 
of the particle with the static chain network. Number of par- 
ticles: 500. Chain volume: 12.5%, 3x3x3 particle (large). 
Lower panel: maximum mean square displacement. Number 
of particles: 500. Chain volume: 12.5%, 3x3x3 particle 
(large). 



